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The  character  I zat i on  of  composites  is  a  continuing  pro¬ 
cess  involving  a  wide  variety  of  materials  with  different 
compositions,  laminate  geometries  and  fabrication  processes. 
This  thesis  examines  bench  test  methods  for  the  determination 
of  the  linear  coefficients  of  thermal  expansion,  residual 
stresses  due  to  fabrication,  and  del  ami  nation  strengths  in 
S2-G 1 ass/Epoxy  and  AS4  Graphite/Epoxy  filament  wound  cylin¬ 
ders.  Results  from  these  tests  are  compared  to  finite  ele¬ 
ment  and  closed  form  solutions.  Test  theory  correlation  is 
generally  good.  The  slopes  of  temperature-strain  curves  are 
closely  approximated  by  the  linear  coefficient  of  thermal 
expansion  calculated  for  the  laminates.  Predicted  values  of 
residual  stress  are  approximately  equal  to  those  determined 
by  a  simple  bench  test.  Measured  del  ami  nation  strengths  were 
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CHAPTER  1 


INTRODUCTION 


1 . 1  Background 

A  composite  material  can  be  defined  as  a  materials  system 
composed  of  a  mixture  or  combination  of  two  or  more  macrocon¬ 
stituents  differing  in  form  and/or  material  composition  and 
that  are  essentially  insoluble  in  each  other.  Generally,  we 
consider  advanced  composites  to  be  fiber-matrix  composites. 
Even  within  this  classification  of  composites  there  exist 
millions  of  possible  final  usable  materials.  Many  complex 
processing  variables  control  the  quality  and  reliability  or 
uniformity  of  composites.  They  relate  both  to  the  quality  of 
the  individual  constituent  elements  and,  often  more  impor¬ 
tant,  to  that  involved  when  the  constituents  are  combined 
into  the  finished  composite.  In  many  fiber  composites,  uni¬ 
formity  in  orientation  of  the  fibers,  control  of  the  matrix 
content,  quality  of  the  inter facial  bonds,  and  many  other 
processing  variables  influence  the  composite's  performance 
d i rect 1 y . 

Even  if  standards  of  quality  and  uniformity  were  perfect, 
there  are  still  a  vast  number  of  possible  fiber  and  matrix 
materials  to  choose  from.  Further,  for  every  fiber-matrix 
pair  there  are  an  infinite  number  of  laminate  geometries. 
Given  all  these  different  materials,  character i zat i on  of  com¬ 
posite  materials  In  general  becomes  a  challenging  task 


without  an  adequate  unifying  theory. 


To  determine  the  strength  and  other  properties  of  a  pat — 
ticular  composite  researchers  are  forced  to  either  compare  it 
to  a  similar  material  for  which  characterization  testing  has 
already  been  accomplished  or  perform  testing  ourselves. 

Bert'  recognized  that  simple  bench  tests  are  desirable  in 
determining  these  properties.  This  thesis  investigates  meth¬ 
ods  for  the  determination  of  the  linear  coefficients  of  ther¬ 
mal  expansion  (Chapter  2).  residual  stresses  due  to  fabrica¬ 
tion  (Chapter  3),  and  delamination  strengths  (Chapter  4)  in 
S2-G 1 ass/Epoxy  and  AS4  Graphite/Epoxy  filament  wound  cylin¬ 
drical  composites. 

The  prediction  of  laminate  properties  from  known  or 
assumed  fiber  and  matrix  properties  is  also  a  difficult  task. 
Another  goal  of  this  paper  is  to  determine  the  applicability 
of  the  proprietary  program  AN5Y5,  a  well-known  general  pur¬ 
pose  finite  element  program,  to  the  problems  covered.  An 
educational  version  of  ANSYS  was  available  at  UNL  on  the  VAX 
mainframe  system.  Finite  element  results  are  then  compared 
to  closed  form  solutions  using  other  locally  available  pro¬ 
grams  . 

1.3  Specimens  and  Properties 

Experimental  work  began  with  tests  being  conducted  in  the 

2 

Engineering  Mechanics  lab  at  UNL  by  Foral  et . a  I  .  with  f i  la¬ 
ment  wound  cylinders.  Four  different  types  of  laminates  were 


tested;  two  compression  and  two  tension.  Al I  tubes  were 


wound  on  a  2.78  inch  diameter  mandrel  with  the  compression 
specimens  being  cut  to  approx (mate  1 y  3"  in  length  and  the 
tension  specimens  being  cut  to  approx imate 1 y  14".  Table  1.1 
(next  page)  shows  the  various  laminates,  thicknesses,  and 
wind  angles.  "Inner"  refers  to  the  all  S2  Glass/Epoxy  lami¬ 
nates  while  "outer"  refers  to  the  laminates  with  graphite 
helicals.  Information  concerning  the  "outer"  tension  tests 
significant  to  this  thesis  is  the  location  of  de 1  ami nat 1  on . 
Large  areas  in  the  test  section  of  the  cylindrical  specimens 
delaminated  at  failure  consistently  between  the  gra¬ 
phite/epoxy  and  glass/epoxy  helical  layers  (Table  1.1). 

Some  of  the  specimens  being  tested  in  compression  were 
cooled  or  heated  to  certain  temperatures  prescribed  for  the 
particular  compression  tests.  As  the  opportunity  presented 
itself,  strain  data  was  recorded  during  the  change  in  temper¬ 
ature  using  an  IBM  PC  Portable  computer.  Measurements  Group 
2100  System  Signal  Conditioner  and  Amplifier,  Metrabyte 
Oasn-8  A/D  Converter,  a  program  cal  led  DASHSOfT  written  by 
Gilbreath3,  and  an  Omega  2 1 68A  Digital  Thermocouple  device. 
This  data  was  then  used  to  determine  a  coefficient  of  thermal 
expansion  and  residual  stresses  in  the  specimens  in  Chapters 
2  and  3 . 

Although  four  specimen  types  were  used,  this  paper  usu¬ 
ally  discusses  only  one  specimen  type  in  detail.  Material 

4 

properties  used  were  taken  from  Chamis  with  properties  for 
♦  /-  theta  sets  of  he  I  i ca I  pi  ies  being  calculated  with  lami- 
nat  e  t  heor  y  using  three  di  mens  i  ona  I  cons  i  derat  i  ons<> .  Table 
1.2  shows  the  properties  used  in  this  thesis. 


INNER  LAMINATES 


1 


TENSION 

PLY 

1 

2 

3 

4 

5 


THICKNESS 

0.01 

0.00746 

0.00765 

0.00765 

0.00746 

0.04022 


ANGLE  FROM  AXIAL 
0 
90 
+20 
-20 
90 


FIBER/MATRIX 
SCRIM  CLOTH 
S2  GLASS/EPOXY 


COMPRESSION 

(PLY  #1  ABOVE  THEN  REPEAT  PLIES  2-6  THREE 
TIMES) 


OUTER  LAMINATES 


1 

TENSION 

0.00746 

COMPRESSION 

0.0224 

90 

S2 

GLASS/EPOXY 

2 

0.00746 

0.0224 

+  20 

AS4  GRAPHITE/EPOXY 

3 

0.00746 

0.0224 

-20 

t» 

4 

0.00765 

0.0229 

+20 

S2 

GLASS/EPOXY 

5 

0.00765 

0.0229 

-20 

ft 

6 

0.00746 

0.0224 

90 

ft 

0.04514 

0. 1354 

Table  1.1  --  Cylindrical  specimen  configurations. 
Thicknesses  in  inches.  Layers  listed  in  order  of 
inside  to  outside  through  thickness. 
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G 1  2 

(psi  ) 

0.602 

2.632 

G 1  3 
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0.602 
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(see  Tab 
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( see  Tab 
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20.15 
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CHAPTER  2  -  COEFFICIENT  OF  THERMAL  EXPANSION 


2.1  Introduction 

The  temperature  related  experimental  results  were  used  to 
determine  a  coefficient  of  thermal  expansion  (CTE)  for  both 
the  inner  and  outer  compression  specimens.  These  results 
were  then  compared  to  results  obtained  in  closed  form  and 
finite  element  solutions  and  finally  to  results  found  in 
recent  literature. 

Key  F i gures : 

Figures  2.5  and  2.6  present  measured  and  corrected  ther¬ 
mal  expansion  test  data  (heavy  solid  line),  and  the  theoreti¬ 
cal  correlation  (dashed  line). 

2.2  Experimental  Results 


2.2.1  Test i ng 

Figures  2.1  and  2.2  show  the  actual  readings  taken  from 
the  test  described  on  page  3  for  the  outer  and  Inner  spe- 
c i mens  respect ively.  In  Figure  2.1,  the  varl ous  symbo is  in 
the  pre-75  degree  range  and  those  in  the  post-75  degree  range 
represent  different  specimens.  The  same  applies  to  Figure 

2.2  for  the  whole  range  of  temperatures,  since  only  five 
inner  specimens  were  tested  in  this  manner  (the  repetition  of 
symbols  is  due  to  a  limitation  of  LOTUS  123  software).  Keep 
in  mind  that  the  nonlinearity  of  these  graphs  should  be 
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decreased  by  the  superposition  of  compensation  gage  data. 


Since  there  were  no  compensation  gages  used  in  the  origi¬ 
nal  testing,  it  was  necessary  to  run  additional  compensation 
tests  in  order  to  determine  the  coefficient  of  thermal  expan¬ 
sion.  Preferably,  one  would  use  identical  strain  gages 
mounted  on  low-expansion  National  Bureau  of  Standards  (NBS) 

standard  reference  material  of  fused  silica  concurrent  with 
S  7  8 

actual  testing.  ’  ’  Also,  it  is  desirable  to  heat  the  speci¬ 
men  in  steps  allowing  thermal  equilibrium  between  successive 
steps.  This  was  not  possible  due  to  a  prescribed  temperature 
rate  for  the  original  tests  of  5°F  per  minute. 

The  gages  used  for  compensation  measurement  were  identi¬ 
cal  to  those  used  on  the  original  specimens  (see  the  Appendix 
for  gage  information).  Aluminum,  a  readily  available  material 
with  a  known  linear  coefficient  of  thermal  expansion,  was 
used.  Two  gages  were  mounted  (one  centered  on  each  side)  on  a 
2"x5"xl/8"  6061  T6  Aluminum  specimen  using  the  mounting  meth¬ 
ods  described  in  Micro-Measurements  Instructional  Bulletins 
B-129-29  and  8-127-9*^.  The  IBM  PC  Portable,  Measurements 
Group  2100  System  Signal  Conditioner  and  Amplifier,  Metrabyte 
Dash-8  A/D  Converter,  DASHSOFT^  program,  and  Omega  2 1 68A 
digital  thermocouple  device  described  in  Chapter  1  were  used. 
The  specimen  was  placed  in  a  Precision  Scientific  Company 
2700  watt  oven  (range  i  from  room  temperature  to  356°F).  The 
temperature  was  increased  at  approx (mate  I y  five  degrees  Fah¬ 
renheit  per  minute  as  witnessed  by  the  time  vs.  temperature 
graph  of  1-  igure  2.3  (approximately  the  same  rate  as  the  or  i  g- 
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Figure  2.3 


TIME  (MIN) 

--  Time  vs.  Temperature  for  Compensation  Test 


6061  T6  ALUMINUM  COMPENSATION  GACE 


for  the  compensation  test  .  Note  how  the  graph  remains 
linear  during  "unloading." 


2.2.2  Theory  and  Reduction  of  Data 

A  linear  regression  analysis**  was  then  performed  on  this 
data  to  arrive  at  a  linear  coefficient  of  expansion  due  to 
both  the  gage  and  the  aluminum  specimen.  To  arrive  at  a  CTE 
for  the  strain  gage  only,  the  methods  described  by  Poore  and 
Kesterson^  and  Measurements  Group  Tech  Note  513^  were  used 
with  slight  modifications  to  allow  for  the  fact  that  compen¬ 
sation  gage  testing  was  done  after  the  original  tests. 

We  have. 


CTE 


compos i te 


where. 


CTE 

measured 

CTE 

measured 


+  CTE 

compensat ion 
+  *CTEaluminum 


CTE*  ) 


CTEmeasured  =  slope  of  stra i n-temperature  graph  for 

m  original  composite  specimens. 

CTE  =  slope  of  stra i n-temperature  graph  for  the 

aluminum  compensation  material. 

CTE  .  =  the  CTE  known  for  6061  T6  aluminum, 

a l um  f  num 


The  slope  of  the  stra i n-temperature  graph  for  the  alumi¬ 
num  alloy  (Figure  2.4)  was  found  to  be  7.5  microstrain  per 


degree  F.  Subtracting  this  from  a  "known"  value  for  the  CTE 

22 


of  the  compensation  material  of  13.1  we  have  CTE 

x  com¬ 
pensation  =  ^*6  •  Tables  2.1  and  2.2  show,  for  the  outer 

specimens,  the  measured  CTE  and  the  corrected  CTE  which  rep¬ 
resents  the  CTE  of  the  composite  cylinder  in  the  axial  and 
hoop  directions  respectively.  Starred  specimens  are  those 


with  a  coefficient  of  correlation  greater  than  0.998  calcu- 


lated  for  CTE 


measured 


Values  of  CTE  for  highly  non-linear 


data  are  measured  on  the  first  linear  portion  as  the  tempera¬ 
ture  drops  below  room  temperature  as  indicated  in  Figure  2-1. 
The  mean  (AVE),  standard  deviation  (SDEV),  and  coefficient  of 
variation  (COV)  are  calculated  using  only  values  from  these 
starred  specimens. 

Referring  back  to  Figure  2.1  for  the  outer  specimens,  the 
high  temperature  data  is  obviously  fairly  linear  and  linear 
coefficients  of  expansion  were  easily  measured.  However,  the 
1 ow-temperature  specimens  show  a  very  non-linear  relationship 
of  strain  with  temperature.  This  non-linearity  makes  the 
determination  of  a  linear  CTE  difficult  as  Table  2.2  shows  a 
high  coefficient  of  variation  among  the  axial  tubes.  (The 
"hoop"  direction  gages  are  included  in  the  tables.  Their 
graphs  are  not  shown  for  brevity,  however.  Although  the 
plots  were  fairly  linear,  little  can  be  concluded  from  one 
specimen  in  each  temperature  range.) 


SPECIMEN 

DIRECTION 

CTE 

measured 

compos i te 

« 

2-6-6 

AXIAL 

-3.96 

1 .64 

» 

2-4-3-3 

AXIAL 

-3.83 

1 .77 

• 

1 -2-2-4 

AXIAL 

-4.14 

1 .46 

« 

1  -  1  -3-3 

AXIAL 

-4.08 

1 .52 

« 

3-2-4-2 

AXIAL 

-3.69 

1 .91 

• 

3-4-3-2 

AXIAL 

-4.16 

1  .  44 

AVE 

=  -3.98 

1.62 

SDEV 

=  0.19 

0.  19 

COV 

=  0.05 

0.12 

* 

2-6-6 

HOOP 

2.13 

7.73 

Table  2.1 
Hot-Dry  Outer 

^/5f. 


Coefficients  of  Thermal  Expansion  for 
Specimens.  All  dimensional  values  are 


i  n 


SPECIMEN 

DIRECTION 

CTE 

measured 

CTE  .  . 

compos i te 

3-4-3-3 

AXIAL 

-4.30 

1 .30 

3-4-4-2 

AXIAL 

-3.36 

2.24 

3-4-5 

AXIAL 

-2.76 

2.84 

3-4-6 

AXIAL 

-2.36 

3.24 

AVE 

=  -3.19 

2.41 

SDEV 

=  0.84 

0.84 

COV 

=  0.26 

0.35 

*  3-4-6 

HOOP 

3.61 

8.66 

Table  2.2  --  Coefficients  of  Thermal  Expansion  for 
Low-Tempergture  Outer  Specimens.  All  dimensional  val¬ 
ues  are F  • 

Tables  2.3  and  2.4  show  this  information  for  the  inner 
specimens  in  h i gh-temperatur e  and  low-temperature  tests 
respectively.  The  high  non-linearity  in  the  inner  specimens 
occurred  in  the  high-temperature  range  while  the  low- 
temperatire  range  results  are  fairly  linear  as  witnessed  by 
Figure  2.2.  As  stated  earlier,  we  expect  the  non-linearity 
of  Figures  2.1  and  2.2  to  decrease  when  we  adjust  them  to 
include  the  effect  of  the  compensation  gage. 


SPECIMEN 


DIRECTION  CTE 


measured 


CTE 


compos i te 


1-2-1-1 

AXIAL 

-5.92 

-0.32 

3-3-4-2 

AXIAL 

-2.63 

2.97 

5-1-2-1 

AXIAL 

-3.46 

2.14 

AVE  = 

-4.00 

1 .60 

SDEV  = 

1.71 

1.71 

COV  = 

0.43 

1 .07 

5-1-2-1 

HOOP 

( INCONCLUSIVE) 

Table  2.3  --  Coefficients  of  Thermal  Expansion  for 
Hot-Dry  Inner  Specimens.  All  dimensional  values  are  in 
^/°F- 


SPECIMEN 


DIRECTION 


CTE 


measured 


CTE 


compos i te 


*  3-3-5 

AXIAL 

6.67 

12.27 

*  3-4-5 

AXIAL 

7.82 

13.42 

AVE  =  7.24 

12.84 

SDE V  =  0.81 

0.81 

COV  =  0.11 

0.06 

•  3-4-5 

HOOP 

7.84 

13.44 

Table  2.4  --  Coefficients  of  Thermal  Expansion  for 
Low-Tempergture  Inner  Specimens.  All  dimensional  val¬ 
ues  are  z*/  r  ■ 


The  equation  for  the  linearized  form  of  Figure  2.4  is 
STRA I N=5 . 6*TEMPERATURE-464 .  If,  instead  of  measuring  the  CTE 
of  the  specimens  before  adding  the  effect  of  the  compensation 
gage,  this  equation  is  superimposed  onto  Figures  2.1  and  2.2; 
the  non-linearity  can  be  reduced  and  the  actual  CTE  of  the 
material  can  be  measured  directly.  Consequently,  Figures  2.5 
and  2.6  are  simply  this  superposition.  They  indicate  the 
actual  strain  of  the  cylinders  under  change  in  temperature. 

The  data  points  still  Indicate  a  significant  non-linearity  despit 
an  obvious  reduction  from  Figures  2.1  and  2.2.  This  non¬ 
linearity  represents  a  problem  in  determining  a  single  linear 
coefficient  of  thermal  expansion  for  the  test  materials.  The 
solid  curved  lines  In  Figures  2.5  and  2.6  Indicate  a  cubic  spline 
analysis  using  five-point  smoothing^.  Both  curves  have 
coefficients  of  correlation  greater  than  0.9990.  The  dashed 
lines  in  Figures  2.5  and  2.6  represent  the  data  predicted  by 
the  closed  form  and  finite  element  solutions  in  Sections  2.3 


and  2.4. 


Finally,  the  heavy  solid  line  indicates  the  pre- 


dieted  thermal  expansion  after  compensation  adjustments  dis¬ 
cussed  in  the  following  paragraphs. 


2.2.3  Compensation  Adjustment 

The  non-linearity  of  both  inners  and  outers  is  consistent 
among  the  data  available.  One  possibility  may  be  that  the 
compensation  measurement  is  not  accurate  enough.  A  decrease 


in  the  compensation  value  (5.6 


}F )  would  decrease  the 


amount  of  reduction  in  non-linearity  while  an  increase  in  the 
compensation  would  decrease  the  non-linearity  of  the  final 


material  strain-temperature  curve. 


Therefore,  if  CTE 


compensat i on 


were,  say,  5.7  >*/  F ,  the 


curve  of  Figures  2.5  and  2.6  would  be  more  linear,  thus 

decreasing  the  change  in  CTE  with  change  in  temperature. 

Since  a  value  of  CTE  for  aluminum  of  13.1  ^///°F  was  used,  an 

error  of  only  0.757.  in  this  value  can  significantly  change 

the  data.  This  is  why  it  is  preferable  to  use  a  compensation 

material  with  a  very  low  coefficient  of  expansion,  such  as 

the  fused  silica  mentioned  earlier.  If  CTE  .  .  is 

compensat i on 

increased  until  Figures  2.5  and  2.6  become  nearly  linear  over 
a  range  not  including  the  extreme  ends  (using  LOTUS  123),  we 


arrive  at  CTE  „  . ,  =5.67  and  a  linear  coefficient  of 

compensat  f  on 

expansion  of  the  test  material  of  9.72 ^></0F  and  2.63^/°F  for 
the  inner  and  outer  specimens  respectively  (see  Table  2.9  on 


page  28  for  compared  results).  This  corresponds  to  a  CTE  for 
6061  T6  Aluminum  of  13.17 ^/°F  compared  to  the  original 
assumed  value  of  13.10 i^/t,/°F  from  Military  Handbook  17^. 

This  value  was  consistent  for  both  inner  and  outer  specimens. 


2.3  Closed  Form  Solution 


16 


2.3.1  TDCYL2  Program 

Closed  form  solutions  were  found  for  Inner  and  outer, 
tension  and  compression  specimens  using  the  FORTRAN  program 
TDCYL2  written  by  N.L.  Newhouse^.  TDCYL2  uses  a  closed  form 
method  to  calculate  the  stresses  and  strains  In  the  three 
principal  directions  at  any  radial  position  through  the 
thickness  of  a  composite  cylinder.  TDCYL2  accepts  material 
properties  only  In  the  meridional,  hoop,  and  radial  direc¬ 
tions  (m,h,r)  as  shown  In  Figure  2.7. 

Y 


Figure  2.7  --  Coordinate  systems. 

Therefore,  In  order  to  use  TDCYL2  In  the  determination  of 
thermal  strains,  coefficients  of  thermal  expansion  must  be 
Input  for  each  layer  In  the  meridional,  hoop,  and  radial 
directions.  As  outlined  In  Table  1.2  the  Chamls^  properties 
for  CTE  were  used,  namely. 


S- GLASS /EPOXY 


AS- 4  GRAPHITE /EPOXY 


— i  -  *1  .  c.  c. 

oti  =  20 .  1  5 
<*^  =  20.  15 


I* 

*« 


«.  =  0  . 06  ^n/°F 

oJ= 18.06  " 

o<2=  18.06  ” 


2.3.2  Determination  of  Helical  Ply  Coefficients 

Now,  for  the  transformat i on  of  CTE  from  the  1,2  coordi¬ 
nate  to  the  m,h  coordinate,  we  have. 


“  — 

— —  — 

<VT 

= 

2  2 

m  n  2mn 

2  2 

n  m  -2mn 

«j4T 

<x2AT 

2  2 

-mn  mn  m  -n 

where  m=cos(fl)  and  n=sin(<J>.  This  reduces  to, 

(2-2) 
(2-3) 
(2-4) 


<V  =  «cos2fl  +  «tsln2« 
m  l  c. 

2  2 
orh  =  ^sin  +  «2cos  9 

<vmh  =  -WjCostfs  in^  +  o^cos^sintf 


As  long  as  we  are  dealing  with  0  and  90  degree  plies  we 

can  use  equations  (2-2)  through  (2-4)  since  (2-4)  becomes 

zero  (TDCYL2  does  not  accept  .  as  input).  However,  since 

mh 

we  have  +/-  20  pi ies  we  must  find  a  way  to  cancel  out  .  . 

mn 

This  can  be  done  by  including  a  +  9p}y  along  with  a  -0ply  as 

a  single  layer  in  an  analysis  using  lamination  theory. 

The  problem,  then,  is  to  find  or  ,  or  ,  or  for  a  +/-  theta 

m  n  r 

laminate  given  o^  ,  ^  and  the  stiffness  properties  of  the 


i  nd i v I dua I  plies. 


For  this  ana  lysis,  let  the  coordinate  sys- 


tem  x,y,z  be  congruent  to  m,h,r.  From  general  laminate 


theory  we  have  equations  (2-5),  (2-6),  and  (2-7);  where  [T] 

Is  the  two-dimensional  transformat  Ion  matrix,  [T]  *  Is  its 

Inverse,  [Q.  is  the  stiffness  matrix  for  the  kth  layer, 
i  J  k 

[  (J~  ].  Is  the  stress  vector  for  the  kth  layer,  and  [£"  ]. 

k  x » y  k 

is  the  strain  vector  for  the  kth  layer. 


(2-5) 


[T] 


-1 


(2-6) 


.-1 


[°-x,yik  -  fn  *CQ,jVTH€x,y]k 


(2-7) 


We 


then  define  as  the  matrix  multiplications 


,-1 


[T]  [Qjj]k[T]  and  equation  (2-7),  with  the  thermal  expansion 


terms  added  in,  becomes. 


<rx 

*12 

6 

£  - 
X 

ot  AT 

X 

°y 

Q)2 

Q22 

2*26 

£  - 

y 

ocyAJ 

^xy 

k 

^16 

^26 

2*66 

k 

£ 

xy 

o<XyAT/2 

(2-8) 


For  the  thickness,  t^,  distance  to  midlayer  from  refer¬ 


ence,  z.  ,  and  the  $  matrix  of  each  layer  we  can  define. 


V) 


(2-9) 


A  1  J  ?,(<7f  j)ktk 


B 


1  J 


^(J)ktk2k 


(2-10) 


where  n  Is  the  number  of  layers  and  Afj  and  B^j 

coupling  midlayer  strain,£T  ,  and  curvature,  k, 

o 

running  load,  N^j. 

[N]  =  [A][£To]  +  [  B  ]  [  k  ]  -  [N]T 
where  [N]^  is  th^  running  thermal  load. 

[N|j]t  .  £ 

k-/  ) 

for  1  ,  j  =  1 . 2 ,6 

We  can  evaluate  the  integral  as  follows: 


are  matrices 
to  the 


(2-11) 


(2-12) 


hk 

\QiJ-\j*Td2  "  QiJ^iJAT(hk-hk-l> 

Jhk-1 

=  Q  M  ATt  (2-13) 

Therefore , 

rN(jJT  -  'i,(0|j>k(’'u)kATtk 

=  £(0,J)k(-1J)ktk  <2-10 

From  equation  (2-11)  we  assume  there  is  no  curvature  so 
[B]  =  0.  Therefore  we  can  say, 

[N]  =  [A][£o]  -  [N]T  (2-15) 

Sine*5  we  are  applying  only  a  change  in  temperature,  we 
can  say  [N]=0  and  it  follows  that, 

[A][£o]  =  [N]T  (2-16) 

Now  premultiply  by  [A]  *, 

t£Q]  =  f A]' 1 fN}T  (2-17) 

Now ,  replacing  [ N ] T  in  equat ion  (2-17)  by  equat ion  (2-14) 


and  writing  in  [£_] 


<2-  IB) 


where  («< .  )  =  (<*)  for  a  theta  degree  ply. 

'  J  &  x .  y 

The  Q  terms  are  easily  computed  and,  hence,  the  A. .  and 
'  •  i  J 

[A]  terms  (a  programmable  calculator  was  used  for  the  3x3 
operation.  A  Gaussian  Elimination'1  procedure  was  used  to 
invert  the  A  mat  rix)  .  Now ,  as  sum ing^T=l  degree  F  and  i nput  - 
ting  the  original  values  for  GTE  in  the  1,2  directions  from 

4 

Chamis  ,  CTE's  for  the  x,y  (m,h)  coordinate  system  can  pe  com 
puted.  The  -xy  or  -mh  terms  drop  out  for  +/-theta  laminates. 
Table  2.5  shows  the  results  for  our  two  materials  in  + / -  20 
degree  configurations.  Since  this  theory  applies  only  to  the 
x,y  (m,h)  coordinates,  the  radial  CTf  I s  assumed  to  be  equal 
to  the  hoop  GTE  as  would  be  expected. 


AS4- Graph  I te /Epoxy  S2-G! ass  'Epoxy 


Of 

1  .  b  x 

1  0'J?/F 

3 

.  4  x 

1  o' 

*/F 

K 

c* 

14.3  x 

1  0  /F 

’  9 

.  r  x 

l  0 

'  /F 

V 

14.3  X 

1  0'b/F 

1  w 

.  n  x 

1  0 

f'/F 

Table  2.5  -  -  C  a  I  cu  I  at  ed  values  for  CTf  for  a  ♦  2  o' 

layer  using  laminate  theory. 

A  t  wo  -  d  i  mens  Iona1  anal/eis  using  aener  a  1  I  ami  i  nafp  their, 

?  i 

on  i nd l  v i dua  '  plies  in  t  he  progr am  WA  T  ^  vie'  ft  t  he  r  pxu  i  r 


of  T ah i e  2 . f . 


F I gure  2 . 8 


CTE  vs.  Ply  Orientation  Angle 


ASA- Graph  1 te /Epoxy 


S2-G I  ass /Epoxy 


-1.57  x  10 _“/F 
15.35  x  1  0 b/F 
15.35  x  10  / F 


3.40  x  10  “/F 
18. 74  x  10  b/F 
18.74  x  10  b/F 


Table  2.6  --  Calculated  values  for  CTE  for  ♦/-  20 
degree  laminate  using  two-dimensional  single  ply 
so  1 ut 1  on . 


Using  the  above  method,  CTE  can  be  plotted  against  dif¬ 
ferent  values  of  theta  for  both  a  single  theta  ply  and  a  +/- 
theta  laminate  as  shown  In  Figure  2.8.  Herakovlch^  discusses 


this  phenomena  with  respect  to  other  material  properties  with 


emphasis  on  negative  through-the-thlckness  Poisson's  ratios. 


2.3.3  TDCYL2  Results 


4 

These  properties  of  CTE  and  the  Chamis  material  propet — 
ties  of  Table  1.2  were  then  entered  into  TDCYL2  with  the 
results  shown  in  Table  2.7  for  tension  and  compression,  and 
inner  and  outer  specimens.  (Radial  CTE ' s  for  graphite  and 
glass  helicals  used  were  19.3 ^u/° F  and  20.8 ^,/°F  due  to 
earlier  predictions  of  CTE). 


TENSION  COMPRESSION 

INNER  OUTER  INNER  OUTER 

°r  8.07  ^./F  2.3E  WF  8.89^,/F  2.4  1  ~/F 

«*  7.26  8.73  7.50  9.16 


Table  2.7  --  T0CYL2  Coefficient  of  Therma I  Expan¬ 
sion  output. 

Only  compression  specimen  CTE  results  can  be  compared 

with  experimental  results  since  only  compression  members  were 

tested  for  CTE.  Referring  back  to  Table  2.2  we  see  that  the 

outer  low-temp  specimens  compare  very  we  1  I  but  in  Table  2.4 

the  closed  form  solution  falls  short  of  the  measured  value. 

Remember,  however,  the  problem  of  non- I i near i t y .  The  values 

for  CTE  after  incrementing  CTE  ,  in  Section  2.2  are 

compensation 

very  close  to  those  predicted  by  TDCYL2  in  Table  2.7  with  a 
coefficients  of  variation  of  6.3%  and  6.2%  for  the  inner  and 
outer  laminates,  respectively  (for  comparison,  see  Table  2.9 


on  page  28 ) . 


2.4  Finite  Element  Solution 


2.4.1  Description  of  AN5YS 

For  a  finite  element  solution,  the  proprietary  program 
ANSYS,  version  4.2,  available  on  the  UNL  VAX  system  was  used. 
ANSYS  is  a  software  engineering  analysis  system  produced  and 
sold  by  Swanson  Analysis  Systems,  Inc.,  Houston,  Pennsylva¬ 
nia.  ANSYS  capabilities  include  structural  analysis  (static 
and  dynamic;  elastic,  plastic,  creep,  and  swelling;  buckling; 
small  and  large  deflection  theory),  heat  transfer  analysis 
( steady- state  and  transient;  conduction,  convection,  and 
radiation;  phase  change),  uncoupled  thermal-stress  analysis 
(temperatures  from  a  heat  transfer  analysis  available  for 
input  into  a  Structural  analysis),  coupled  thermal  stress, 
and  electromagnetic  field  analysis.  ANSYS  was  originally 
designed  (1970)  from  the  ground  up  mainly  for  isotropic  mate¬ 
rials.  Capabilities  to  handle  anisotropic  and  non-linear 
materials  were  added  later.  The  use  of  ANSYS  in  the  investi¬ 
gation  of  composites  appears  to  be  a  rarely  attempted 
endeavor  and  its  inclusion  in  this  paper  serves  not  only  as  a 
finite  element  tool  to  compare  with  other  solutions  but  also 
as  a  study  of  this  general  purpose  program's  app 1 i cab i 1 i ty  to 
research  in  composite  materials  characterization. 

2.4.2  ANSYS  Input 

Figure  2.9  is  the  ANSYS  input  data  for  the  thermal  analy¬ 
sis  of  the  outer-tens i on  specimens  (Data  for  the  other  speci¬ 
men  configurations  is  not  included  for  brevity).  The  "NL" 


LAYERS 


commands  give  the  material  properties  In  the  form  of  a 


6x6  stiffness  matrix  equivalent  to  the  Q  matrix  discussed 
earlier.  The  "N"  and  "E"  commands  define  the  nodes  and  ele¬ 
ments  respect  1 ve 1 y .  Figure  2.10  shows  the  element  configura 
tion  used  to  model  the  composite  cylinder  for  thermal  analy¬ 
sis.  Scale  is  exaggerated  to  show  how  the  elements  define  a 
cylindrical  shell.  In  the  mode),  each  layer  of  composite  is 
characterized  by  three  to  five  elements  in  a 


/PREP7//T I TLE , TENS  I  ON , OUTER , I NTERLAM I NAR .OUT- 
TEMP.  DAT, TEMPERATURE  RUN 
ET, 1 ,64,0,0,0,0,0,0 

NL , 1 , 1 , 2 . 408E6/NL ,1,7,0. 968E6/NL ,1,13,1. 002E6/NL ,1,37,8- 

. 190E6  NL , 1 , 43 , 0 . 968E6/NL , 1 ,67 , 2 . 408E6/NL , 1 ,91 .0.703E6 

NL, 1, 1G9.0.703E6/NL, 1 , 121 ,0.703E6/NL,6, 1 ,2.408E6 

NL , 6 , 7 , 0 . 998E6/NL , 6 , 13,0.972E6/NL,6,37,2.480E6 

NL,6,43, 1 .573E6/NL, 6,67, 6.91 0E6/NL, 6,91 ,0.692E6 

NL , 6 , 1 09, 1 . 869E6/NL ,6,121,1. 504E6/NL ,7,1,1. 356E6 

NL , 7, 7 ,0. 468E6/NL ,7,13,0.51E6/NL,7,37, 1 .698E6/NL,7,43,2- 

NL,7,67,16.665E6/NL,7,91 .0.472E6/NL, 7, 109.2.632E6 

NL , 7, 121.7.576E6/ALPX, 1,20. 15E-6/ALPY, 1 ,4.22E-6 

ALPZ, 1 ,20. 15E-6/ALPX,6,20.80E-6/ALPY,6, 19.0  IE- 6 

ALPZ, 6 , 3 . 36E-6/ALPX ,7,19. 356E-6/ALPY, 7,14. 254E-6 

ALPZ, 7,-1 .621E-6/CSYS, 1 

N, 1 , 1 .390/N,5, 1 .392487/N.9, 1.394974 

N, 13, 1.39746/N, 17, 1.399947/N.21, 1 .402434/N.25, 1.40492 

N, 29, 1 . 40 7404 /N, 33, 1 . 409894/N, 37 , 1 .41 238/N ,41 , 1 .41493 

N,45, 1 .41748/^,49, 1 .42003/N.53, 1 .42258/N.57, 1.42513 

N.61 , 1 .42768/N.65, 1 . 430 167/N ,69, 1 .432654/N, 73, 1 .43514 

NGEN, 2, 1, 1 ,73,4,0,0, .002/NGEN, 2, 2, 1,74,1,0,0. 125,0 

NROTAT ,1,76, 1/MAT , 1/E, 1 ,2, 4, 3, 5, 6, 8, 7 

E, 5, 6, 8, 7,9, 10, 12, 1 1/E, 9, 10, 12, 1 1, 13, 14, 16, 15/MAT, 7 

E, 13, 14, 16, 15, 17,18,20, 19/E, 17, 18,20, 19,21,22,24,23 

E, 2 1,22, 24, 23, 25, 26, 28, 27/E, 25, 26, 28, 27, 29, 30, 32, 31 

E,29,30, 32,31 ,33,34, 36, 35/E, 33, 34,36, 35,37, 38,40,39/MAT- 

,6  E, 37,38,40,39,41 ,42, 44,43/E, 41 ,42,44,43,45,46,48,47 

E, 45, 46,48,47,49, 50, 52, 5 1/E,49, 50,52,51 ,53,54,56,55 

E.53,54, 56,55,57,58, 60, 59/E, 57, 58, 60, 59,61 ,62,64, 63/MAT- 

, 1  E,61 ,62, 64,63,65, 66, 68, 67/E, 65, 66,68,67, 69, 70, 72, 71 

E, 69,70, 72, 71 ,73, 74, 76, 75/CPS  I ZE,40  CP, 1 ,UX, 1 ,3/CP, 2, UX 

,2,4/CP,3.UX,73,75/CP,4,UX,74,76/CPNGEN,5,UZ,2,76,2 

CPNGEN,6,UY, 1 ,73,4/CPNGEN,7,UY,2,74,4/ITER,-5,5 

D,3,UY,0.,,75,4/0,4,UY,0. ,,76,4/D, 1 ,UZ, 0. , , 75 , 2/TREF , 1 0  - 

0  KTEMP, 0/TUN I F , 0/AFWR I TE/F 1 N I SH// 1 NPUT, 2 7/F I N I  SH 


Figure  2.10  --  ANSYS  input  data  for  thermal  analysis. 


through-the-th i ckness  manner.  This  "stack"  of  layers  Is 
then  rotated  about  the  z,  or  axial,  axis  to  form  a  ring 
which,  with  the  proper  boundary  conditions,  represents  the 
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"test  section"  of  a  test  cylinder.  This  model  works  well  in 


representing  a  long  cylinder  under  thermal  or  tensile  loads 
(A  buckling  analysis  requires  a  model  which  includes  the 
entire  length  of  the  cylinder  which  loses  the  three- 
dimensional  capabilities  of  the  present  model.).  The  "CP" 
and  "D"  commands  define  the  coupled  degrees  of  freedom  and 
the  boundary  conditions  respectively.  Finally,  "TREF"  and 
"TUNIF"  define  the  initial  reference  temperature  and  the 
applied  uniform  temperature. 

Since  the  material  properties  input  is  in  the  form  of  the 
Q  matrix,  this  must  be  calculated  separately.  Coefficients 
of  thermal  expansion  for  the  individual  plies  must  also  be 
calculated  and  input  for  the  axial,  hoop,  and  radial  direc¬ 
tions  thus  requiring  the  same  calculations  used  in  TDCYL2 . 

The  stress  results  from  TDCYL2  and  ANSYS  will  be  dis¬ 
cussed  in  Chapter  3,  Residual  Stresses.  CTE  results  from 


TENSION 

COMPRESSION 

INNER 

OUTER 

INNER 

OUTER 

axial 

hoop 

8.074  ^</F 
7.260 

2.376  ^F 
8.732 

8.887  u/ F 
7.503 

2.408  z//F 
9.156 

Table  2.8  --  ANSYS  Coefficient  of  Thermal  Expansion 
output . 

ANSYS  are  shown  In  Table  2.8.  The  results  are  essentially 
the  same  as  for  the  three-dimensional  closed  form  solution 


using  TDCYL2  (Table  2.7).  The  similarity  of  these  results 


provides  some  Insight  Into  the  app 1 i cabi 1 i ty  of  ANSYS  to  com¬ 
posite  materials  and,  specifically  to  thermal  modeling  of 
cylindrical  laminates. 

2.5  Literature  Results 

The  comparison  of  exact  results  for  our  particular  lami¬ 
nates  is  certainly  not  possible  due  to  the  uniqueness  of  the 
layering  variables.  This  is  often,  if  not  always,  a  problem 
in  the  study  of  applicable  composite  laminates.  The  variety 
of  possible  fiber  and  matrix  materials  and  their  countless 
configurations  is  so  great  that  the  classical  comparison  of 
empirical  data  used  in  the  study  of  metals  is  an  impossible 
task.  We  often  rely  on  trends,  general  observations,  or  pos¬ 
sibly  models  developed  from  data  resulting  from  various  con¬ 
figurations  of  like  materials. 

Freeman  and  Campbell15  studied  the  thermal  expansion  of 
various  graphite  reinforced  composites  under  temperatures 
ranging  from  ~300°F  to  500°F.  They  obtained  results  similar 
to  Figure  2.8  experimentally  thus  indicating  the  variability 
of  CTE  with  various  ply  orientations.  They  also  reported 
CTE's  at  three  different  ranges  of  temperature;  -320  to  75°F, 
75  to  350°F,  and  -320  to  350°F.  CTE's  sometimes  tripled  from 
the  first  to  the  second  range  for  particular  graphite  compos¬ 
ites.  This  change  is  not  as  dramatic  for  graphite  composites 
considering  their  small  CTE  but  could  be  a  more  important 
consideration  for  laminates  with  higher  CTE's  (ie.  tripling  a 
CTE=0.7  compared  to  tripling  a  CTE=4  yu/°F).  Freeman  and 
Campbell  reported  an  increase  in  CTE  with  Increases  in  tern- 


perature . 

Daniel  also  reports  significant  increases  in  thermal 
expansion  coefficients  with  increases  in  temperature  for  gra¬ 
phite/epoxy  and  S-g I  ass/epoxy  unidirectional  composites. 
However,  for  some  laminates  with  both  0  and  +/-  45  plies  the 
CTE  actually  decreased  with  increasing  temperatures .  Suffi¬ 
cient  data  was  not  available  to  make  any  general  conclusions 
about  laminate  conf i gurat i on  and  change  in  CTE  with  tempera¬ 
ture.  This  appears  to  be  a  possible  subject  for  further 
i nvest i gat i on . 

2.6  Conclusions 

Table  2.9  indicates  test  and  theory  coefficients  of 
thermal  expansion  as  presented  in  this  chapter.  A  conclusion 
is  that  linear  coefficients  of  thermal  expansion  for  filament 
wound  laminates  can  be  predicted  using  classical  laminate 
theory . 


LINEAR  COEFFICIENTS  OF  THERMAL  EXPANSION 
for  COMPRESSION  SPECIMENS  (m i crostra 1 n/°F ) 


TEST 

THEORY 

'TER 

INNER 

OUTER 

(INSUFFICIENT  DATA) 


Table  2.9  --  Comparison  of  Test  vs.  Theory  for  Linear 
Coefficients  of  Thermal  Expansion  for  Compression  Spec 


CHAPTER  3  -  RESIDUAL  STRESSES 


3.1  Background 

Residual  stresses  are  a  function  of  many  parameters,  such 
as  type  of  fiber  and  matrix,  fiber  volume  ratio,  ply  orienta¬ 
tion,  curing  temperature  and  other  variables.  They  can  reach 
values  comparable  to  the  transverse  strength  of  the  ply  and 
thus  induce  cracking  of  that  ply  within  the  laminate.  They 
reach  equilibrium  with  interlaminar  shear  stresses  trans¬ 
mitted  from  adjacent  plies  and  can  thus  cause  de I  am i nat i on . 

In  the  design  and  evaluation  of  composite  structures  one  must 
take  residual  stresses  into  account  and  superimpose  them  onto 
those  produced  by  subsequent  externa?  loading. 

Residual  stresses  are  often  called  "fabrication  stresses" 
since  they  result  from  the  fabrication  process.  An  important 
factor  of  this  process  is  the  cure  cycle.  From  the  cure 
cycle  we  can  determine  a  change  in  temperature  from  maximum 
cure  temperature  to  test  temperature.  For  both  inner  and 
outer  specimens,  the  maximum  cure  temperature  was  2 75  +/-10°F 
and  the  typical  test  temperature  was  72  +/-  5°F.  This  yields 
a  change  in  temperature  of  188  to  218°F.  Using  this  change 
in  temperature  we  can  approximate  the  inherent  residual 
stresses  at  the  time  of  testing.  These  residual  stresses  can 
then  be  superimposed  upon  the  actual  test  results  to  yield  a 
more  accurate  picture  of  the  actual  stresses  at  failure. 


As  in  Chapter  2,  experimental  results  are  compared  to 


closed  form  and  finite  element  solutions.  Results  from  the 


TDCYL2  and  ANSYS  programs  will  be  discussed  first,  followed 
by  a  simple  test  procedure  to  determine  residual  stresses  in 
cylindrical  composites. 

Key  F i gure : 

Figure  3.1.b  shows  the  hoop  residual  stresses  as  predicted 
by  TDCYL2  and  ANSYS  as  the  solid  line.  The  dashed  line 
represents  those  residual  stresses  predicted  by  the  bench 
test  discussed  in  Section  3.3. 

3.2  Closed  Form  and  Finite  Element  Solutions 

These  solutions  are  discussed  simultaneously  since  the 
methods  were  covered  in  Chapter  2  and  the  results  from  each 
solution  are  very  similar.  Using  an  average  AT  of  203°F 
+/- 1 5°F ,  Tables  3.1  and  3.2  are  the  results  obtained  from 
TDCYL2  and  ANSYS  respectively  for  the  outermost  fibers  only. 
Using  only  temperature  effects,  these  results  are  our  best 
analytical  prediction  of  residual  stresses.  Since  the 
stresses  in  Tables  3.1  and  3.2  are  for  the  rad i a  1 1 y-outermost 
fibers,  the  radial  stress  is  a  stress  normal  to  the  surface 
which  is  equal  to  zero.  Although  residual  stress  calcula¬ 
tions  were  performed  on  all  four  types  of  specimens,  only 
complete  results  for  the  outer  tension  specimens  will  be 
shown  for  brevity. 

The  three  graphs  of  Figure  3.1  show  the  distribution  of 
radial,  hoop,  and  axial  stresses  through  the  thickness  of  the 
specimen.  Tensile  radial  stresses,  which  would  add  to  a 
delaminating  critical  stress,  are  highest  between  the  glass 
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STRESS 

TENSION 

COMPRESSION 

( ps  i  ) 

INNER  OUTER 

INNER  OUTER 

Table  3.1  --  TDCYL2  output  for  T=-203  F 
va 1 ues  are  in  ps i . 


STRESS 

TENSION 

COMPRESSION 

(PS  i  ) 

INNER  OUTER 

INNER  OUTER 

Table  3.2  --  ANSYS  output  forAT=-203  F.  All  val¬ 
ues  are  i n  ps i . 


helicals  and  the  glass  hoops,  yet  are  very  low  compared 

to  failure  allowables  discussed  in  Chapter  4.  The  hoop  and 

axial  residual  stresses  are  more  significant.  The  only  ten 

sile  hoop  stress  is  4.65  ksi  in  the  glass  helical  layer.  A 

value  of  4.37  ksi  of  this  stress  acts  in  a  transverse  direc 

tion  to  the  glass  fibers.  This  contributes  considerably  to 

23 

an  assumed  failure  stress  of  7  ksi  for  S-g 1 ass/Epoxy  in 
transverse  tension.  More  significantly,  the  axial  residual 
stress  in  the  glass  hoop  layers  is  as  high  as  6.69  ksi  in 
tension.  This  is  very  near  a  failure  stress  and  we  might 
expect  some  matrix  cracking  in  the  specimens  after  cure  or 
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Figure  3.1  —  ANSYS  Residual  Stresses  for  Outer  Tension  Spec. 
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under  small  load  conditions. 


Since  there  appeared  to  be  no  cracking  of  the  matrix  in 
any  of  the  specimens  prior  to  testing,  these  residual 
stresses  may  be  a  high  prediction  or  the  allowable  may  be 
low.  This  shows  the  need,  not  only  for  testing  of  al low- 
ables,  but  also  for  tests  to  measure  residual  stresses.  The 
following  section  will  investigate  a  method  to  measure  resid¬ 
ual  stresses  in  composite  tubes.  Testing  of  allowables 
will  be  discussed  in  Chapter  4. 

3 . 3  Residual  Stre ss  Bench  Test 

26) 

Voyiadjis  et . a  I  .  measured  residual  stresses  in  compos¬ 
ite  cyl inders  by  removing  material  from  the  inner  and  outer 
surfaces  and  measuring  the  resulting  strain  on  the  opposite 
surface.  This  same  idea  can  be  used  more  simply  by  relieving 
the  stress  and  measur i ng  a  change  in  radius  instead  of  a 
change  in  length.  By  measuring  the  change  in  inner  radius 
upon  relief  of  the  residual  stresses  in  the  hoop  direction, 
these  stresses  can  be  calculated.  The  only  unknown  is  the 
change  in  inner  radius. 

To  produce  more  specimens  and  make  cutting  of  specimens 
easier,  3/4",  I",  and  1  1/4"  rings  were  cut  from  an  outer 

tension  cylinder  which  had  not  been  previously  tested.  The 
inner  surface  was  used  for  measurements  because  of  its 
smooth  surface  due  to  its  fabrication  on  a  mandrel  1.  The 
inside  diameter  was  measured  using  inside  calipers  which  were 
then  measured  using  an  outside  micrometer.  This  being 


recorded,  the  ring  specimens  were  cut  once  at  90  degrees 


off  the  recorded  diameter.  This  cut  relieved  the  hoop 
residue:!  stress  in  the  specimen  and  another  measurement  of 
inner  diameter  was  taken  at  the  same  points  of  measurement  as 
before . 

For  six  tests,  two  for  each  width  of  specimen,  the  aver 
age  initial  radius,  R  ,  and  average  final  radius,  R^ ,  were 
1.390  +/-0.002  inches  and  1.3680  +/-  .002  Inches  respec¬ 

tively.  This  corresponds  to  a  change  in  radius  of  -0.022  +  /  - 
0.002  inches.  The  coefficients  of  variation  were  0.367,  and 
0.687.  respectively.  There  was  no  relation  between  width  and 
change  in  radius.  (The  indication  of  error  in  the  readings 
is  necessary  since  a  relatively  sma II  error  results  in  fairly 
large  changes  in  final  stresses.) 

The  change  in  curvature  can  be  computed  as  fol lows, 

A*  =  ( I /R  -  1/Rf)  =  1/1.390  -1/1.368  =  -0.0116  +  /-  0.0021 

O  r 

2  1 

This  Change  in  curvature  can  then  be  input  to  MATL.  to  get 
the  moment ,  M 

test 

M  =  -0.526  +  /-  0.095  in-lb/in 

test 

Now  this  is  compared  to  the  moment  calculated  from  the 
ANSYS  and  TDCYL2  results  for  hoop  stress  in  Figure  3.1.b. 


"predicted  =  °'033 


The  measured  residual  stresses  can  be  calculated  by  multi¬ 
plying  the  values  in  Figure  3.1  by  0.526/0.438  =  1.20.  In 
Figure  3.1.b,  the  heavy  solid  line  is  the  values  predicted  by 
ANSYS,  the  thinner  sol  id  I  i nes  are  the  error  due  to  error  in 
cure  temperature,  the  heavy  dashed  I i ne  is  the  values  pre¬ 
dicted  by  the  bench  test,  and  the  thinner  dashed  I i nes  are 


the  error  in  that  measurement  due  to  micrometer  reading 


3 . 4  Cone  1  us i ons 

The  overlapping  of  the  ranges  in  Figure  3.1  indicates  a 
good  ab i 1 ity  to  predict  residual  stresses  using  TDCYL2  and 
ANSYS.  One  shortcoming  is  that  both  computer  solutions  and 
measurements  of  the  change  in  radius  rely  on  accurate  deter¬ 
mination  of  the  coefficients  of  thermal  expansion  and  other 
material  properties. 

It  has  been  shown  that  radial  residual  stresses  do  not 
add  to  the  tensile  interlaminar  normal  stresses  in  the  region 
where  delamination  occurred  in  the  original  test  specimens 
(see  p.  3),  but  actually  decrease  them  slightly. 

Fairly  significant  residual  stresses  have  been  shown  in 
the  transverse  direction  in  glass  hoop  and  helical  plies, 
possibly  resulting  in  matrix  cracking  during  cool  down  from 
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E  ora  I  and  G i  1  dr  eat  h ' s  cont  i nuing  r  e  sear  ch  i n  modes  of 
failure  in  composite  materials  indicates,  along  with  current 
literature  ,  a  growing  concern  about  delamination 

failures.  The  modes  of  failure  in  filament  wound  composites 
are  particularly  complex.  Both  predicting  the  type  and 
sequence  of  f a i lures  and  analyzing  the  actual  fai lures  which 
tooi<  place  in  a  closely  observed  loading  are  difficult  yet 
necessary  for  the  character izat ion  of  the  materia1. 

The  purpose  of  this  chapter  is  to  investigate  finite  ele¬ 
ment  solutions  compared  to  closed  form  solutions  of  the 
three-dimensional  state  of  stress  in  a  filament  wound  compos¬ 
ite  cylinder  under  axial  tension.  Those  interlaminar  stresses 
contributing  specifically  to  delamination  failures  will  be 
discussed.  A  discussion  of  the  problem  of  discontinuities  at 
the  ends  of  the  cylindrical  test  specimen  caused  by  static 
test  grips  in  relation  to  delamination  causing  stresses  is 
included.  Finally,  experimental  determination  of  the  allow¬ 
ables  which  govern  the  stresses  at  which  del  ami  nation  fail- 
ures  occur  wi I l  be  attempted. 

In  his  programs  TUBE  and  TDCYP,  Gilbreatn'  analyzes  the 
thr ee-d i mens i ona I  state  of  stress  in  the  test  section  of  the 
tubular  specimen  (that  section  of  the  cylinder  not  affected 
by  the  discontinuities  at  the  ends).  Gilbreath  uses  the  "exact" 


method  developed  by  Waltz  and  Vinson  .  The  method  involves 
five  differential  equations  for  each  distinct  layer  in  a  com 
posite  cylinder.  Gilbreath's  solution  is  compared  to  the 
finite  element  solutions  discussed  in  this  paper. 


Measurement  of  Delamination  Stresses 


4.2.1  ANSYS  Input 

Again,  ANSYS  was  used  for  the  finite  element  solution. 
The  physical  model  is  identical  to  that  used  for  the  deter¬ 
mination  of  the  CTE  and  residual  stresses  as  in  Figure  2.10. 
Input  data  is  similar  to  Figure  2.9.  No  temperature  loads 
were  imposed  and  a  constant  positive  unit  axial  strain  was 
imposed  through  the  "D"  boundary  condition  commands. 


4.2.2  ANSYS  Results 

Figures  4.1  and  4.2  are  the  Inplane  and  out-of-plane 
stress  states  through  the  thickness  of  the  outer  tension 
laminate  for  an  average  tested  stra'n  of  1.57..  In  all 
graphs,  the  squares  and  pluses  represent  ANSYS  results  while 
diamonds  and  triangles  represent  the  closed  form  TDCYL2 
results  obtained  from  the  residual  stress  analysis  discussed 
in  Chapter  3.  The  abbreviations  "glO",  "gl20",  and  "gr20" 
represent  the  glass  hoops,  glass  helicals,  and  graphite  heli¬ 
cals  respectively.  The  vertical  solid  lines  in  the  graphs 
represent  the  interlayer  boundaries  between  these  three  dif¬ 
ferent  types  of  windings  and  materials.  The  dashed  lines  are 
the  boundaries  between  + 6  and  -6  helical  plies. 


R&cAiuS 


Figure  4.1  —  Inplane  stress  results  from  ANSYS 
(squares  and  plusses)  and  TDCYL2  (diamonds  and 
triangles) . 


Considering  only  ANSYS  results  at  this  time,  the  squares 
In  the  top  graph  of  Figure  4.1  Indicate  the  hoop,  or  circum¬ 
ferential,  stresses  as  predicted  by  ANSYS  at  the  time  of  mas¬ 
sive  failure  In  the  specimen  (as  opposed  to  matrix  cracking 
which  occurs  somewhat  earlier).  The  plusses,  then,  represent 


Ffgure  4.2  —  Out-of-Plane  stress  results  from 
ANSYS  (solid  lines)  and  TDCYL2  (dashed  lines). 

the  meridional,  or  axial,  stresses  at  failure.  As 
expected,  these  stresses  are  highest  in  the  graphite  helicals 
and  lowest  in  the  glass  hoops.  The  lower  graph  of  Figure  4.1 
is  the  inplane  shear  ( cl rcumferent ia 1 -ax ial )  through  the 
thickness.  The  actual  analysts  of  these  stresses  with  the 
failure  allowables  available  in  literature  and  the  prediction 
of  modes  of  failure  is  left  to  Fora  1  and  Gilbreath  in  their 
continuing  research  and  upcoming  publications. 

4.2.3  Radial  Del  ami  nation  Problem 

The  interest  here  is  in  the  out-of-plane  radial  stresses 
as  predicted  by  ANSYS  in  the  graph  of  Figure  4.2.  In  the 
experimental  work,  del  ami  nation  failures  were  observed  to 
occur  consistently  in  the  interlaminar  region  between  the 
graphtte  helicals  and  glass  helicals.  The  purpose  here  Is  to 
characterize  the  radial  stresses  present  at  this  point  and 
determine  their  ability  to  cause  failure  in  tension  normal  to 


the  surfaces  of  the  plies. 

To  obtain  the  best  prediction  of  the  actual  stress  state 
we  need  to  superimpose  the  residual  stresses  due  to  fabrica¬ 
tion  determined  in  Chapter  3.  Figure  4.3  is  this  superposi¬ 
tion  for  radial  stress  through  the  thickness. 

From  Figure  4.3,  the  radial  stress  between  the  glass 
helicals  and  the  outer  glass  hoop  ply  is  larger  than  between 


the  glass  and  graphite  helicals  where  failure  occurred.  This 
may  be  evidence  that  a  pure  radial  stress  could  not  have 
caused  the  del  ami  nation  failure. 


RADIUS  ( i nches ) 

Figure  4.3  --  Radial  stresses  from  ANSYS  and  TDCYL2 
with  residual  stressses  included. 

2 

This  problem  leads  to  Foral  and  Gilbreath's  work  in 
failure  mode  determination.  There  may  necessarily  be  a  more 
complex  mode  of  failure  Involving  many  more  variables. 


>’.*•  V-  .'-V  \ 


Fora  I  s  work  in  "pul  lout"  of  fibers  may  give  insight  into 
this  problem. 


4.2.4  End  Discontinuities 

27 

According  to  Waltz  and  Vinson  ,  very  little  change  is 


noted  in  the  interlaminar  norma  I  stress  (actually  a  slight 
decrease)  in  moving  from  the  test  section  of  the  cylinder  to 
the  end  discontinuities  (loading  grips). 


The  failures  observed  in  the  cylinders  discussed  in  this 
paper  do  not  appear  to  have  progressed  into  the  region  of 
load  or  structural  discontinuity.  Using  a  relation  by 


Vinson  and  Slerakowski  we  can  approximate  the  distance,  d, 
from  the  grip  in  the  specimen.  The  relation  for  d  is. 


d  =  4  7  Rh(Dx/D0)  ' c  (4-1 ) 

where  R  is  the  average  radius  of  the  cylinder,  h  (s  the 
thickness,  is  a  term  for  the  axial  stiffness,  and  D  is  a 

term  for  the  circumferential  stiffness  (eqn.s  4-2  and  4-3). 


Dx  =£Qx  ^k/12  +  *k^ 


°©  -2^  <tk3/>2  ♦  ttdl2 


(4-2) 


(4-3) 


For  the  outer  tension  specimens  being  discussed  here,  d 


is  evaluated  to  be  0.92  Inches.  An  attempt  was  made  to  model 
this  region  in  ANSYS.  To  maintain  the  appropriate  boundary 
conditions  it  Is  necessary  to  model  the  portion  of  the  cylin¬ 


der  from  the  end  boundary  conditions  to  a  point  in  the  mem¬ 
brane  test  section.  Ideally,  the  entire  ha  1 f  of  the  cylinder 


...  »  a  »  _  a  _  a  _  .*  _ 


would  be  modeled.  However,  just  one  and  one-half  inches  of 
cyl inder  length  resulted  in  more  elements  and  nodes  than  the 
university  version  of  ANSYS  could  maintain.  Very  small  ele¬ 
ments  are  necessary  to  model  the  various  layers  and  an  aspect 
ratio  as  high  as  1:10  was  attempted  with  no  success.  Sub¬ 
structuring  methods  were  also  done  in  an  attempt  to  reduce 
the  number  of  elements.  Some  ground  was  gained  in  this 
attempt.  Computer  time  and  restraints  make  work  with  this 
program  very  tedious  and  expensive. 

Another  method  of  finding  these  out-of-plane  shear 

stresses  would  be  to  solve  the  5*6=30  simultaneous  differen- 
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tial  equations  (Waltz-Vinson  )  in  full  form.  Programs  such 
as  TUTSIM  and  TK  SOLVER  could  be  used  to  solve  these  equa¬ 
tions.  This  work  was  out  of  the  realm  of  this  paper  but  may 
prove  to  be  worthwhile  in  continued  research. 

Regardless  of  any  future  determination  of  the  mode  of 
failure  that  caused  the  delamination,  characterizing  the 
radial  strength  of  these  particular  filament  wound  specimens 
is  necessary.  This  is  the  purpose  of  the  next  section. 


4.3  Allowables  Testing 

The  need  for  bench  tests  to  determine  out-of-plane  normal 
strengths  of  composite  laminate  cylinders  has  been  estab¬ 
lished.  Two  tests  were  developed  and  run  in  a  step  toward 
this  goal.  (Refer  to  Figures  4.6  and  4.15  for  diagrams  of 
the  apparatus  for  each  test.) 


4.3.1  Half-Ring  Bench  Test 


The  following  is  the  process  and  the  calculations  used  to 
determine  the  radial  stresses  given  a  measured  hoop  strain 
due  to  a  change  in  curvature  in  the  hoop  direction.  The 
physical  process  of  measuring  this  strain  is  described  in  the 
next  section. 


4 . 3 . 1 . 1  Theory 

The  rad i a  1  stress  will  be  determined  for  the  interlaminar 

region  between  the  graphite  and  glass  he  I  I ca 1 s  since  this  is 

where  f a i lure  was  observed  to  occur.  Using  the  program 

2  1 

MATL  .  a  unit  moment  was  input  to  calculate,  using  general 

l am i nat e  theory ,  the  inside  strain  ano  change  in  curvature  on 

an  outer  tension  cylindrical  specimen. 

(Note:  Assume  al l  values  have  a  "y"  subscript  to  denote 

the  hoop  direction.  The  subscripts  "act”  and  "mat'”  refer  to 

the  actual  measured  values  and  the  values  calculated  by  the 
2  I 

program  MAIL  respect i ve  '  y  .  ) 

!  NPt  J  T:  M  =-!  in-  lb'  m 

mat  i 

OUTPUT:  { 

C  mat  I 

dr 

mat  I 

By  the  I  i near  re ' at i on?r i p  bet  *een  app I  i ed  moment  and 
resulting  strain. 


Mact^Mmat I  ~  ^ ar*  ^ mat  I 
o'ving  for  the  actual  moment,  M 


act 


(  4  -  4  i 


By  the 


Mact  ’  Sct^mat  I 
I  i near  relationship  between  moment 

4K  MK  M  .  /M  i 


(4  S  ) 

and  cur  vat  ur  e , 

(  4  -  f  ) 


and  solving  for  actual  change  in  curvature, 


act  ’ 


4K  ^  =  -AK  .  M  ^ 
act  mat  1  act 


4K  ,=4K  ,  ,  *  ^  ./£  .  . 

act  mat  I  act  cmatl 


(4-7) 

Now  substituting  (4-5)  into  (4-7), 

(4-8) 

The  actual  curvature  can  now  be  used  to  determine  the 
radius  of  the  specimen  at  any  measured  value  of  strain  at  the 
point  of  placement  of  the  strain  gage.  From  the  definition 
of  curvature, 

(4-9) 


4Kact  ’  ,/Ro  -  '"V 


where  Rq  and  R^.  are  the  initial  and  final  inner  radii 


respectively.  Solving  for  R^, 


Rf  -  V(1-AKactRo> 


(4-10) 

The  radius  of  curvature  is  now  known  and  only  the  hoop 
stresses  are  needed  to  find  the  radial  stress  as  shown  in 

2  I 

Figure  4.4  .  From  the  strain  distribution  taken  from  MATL 
we  can  get  the  average  strains  for  the  two  layers  involved  as 
percentages  of  the  measured  hoop  strains.  Let  a  and  b  be 
these  percentages  in  decimal  form  for  the  glass  and  graphite 
layers  respectively.  Thus, 


^AVEgl  =  a€act 


^AVEgr  "  b*act 


(4-11) 

where  "gl"  and  "gr"  refer  to  the  inside  glass  hoop  ply 
and  the  graphite  helical  plies  respectively.  From  general 
laminate  theory,  the  hoop  stresses  in  the  layers  can  now  be 
so  I ved  for , 


^gl  =  ^22g1^AVEgl 


°gr  =  Q22gr^AVEgr 


(4-12) 


where  Q22g,  and  Q22g(_ 


are  constant  material  values  taken 


layers  and, 

RAVE  ■  Of  *  <v  *  V>'2  <“-|7> 

This  radial  stress  can  then  be  calculated  once  the  strain  at 
failure  Is  measured. 

4.3. 1.2  Test  Apparatus 

Figure  4.6  Is  a  schematic  of  the  testing  Jig  and  appara¬ 
tus  used  to  produce  a  y-moment  (and  a  transverse  shear  as 
discussed  later)  and,  thus,  a  strain  on  the  inner  layer  of 


Figure  4.6  --  Test  Jig  for  Half-Ring  Delaminatlon  Test. 

the  half-ring  specimen.  Grip  block  "A"  Is  securely  fas¬ 
tened  to  plate  "0"  which  Is  free  to  translate  and  rotate 
about  Its  plane.  Block  "B"  is  free  to  translate  away  from 
block  "A"  as  the  specimen  flattens  out  under  the  load  at 
point  "E".  Block  "C”  was  used  simply  as  a  spacer  and  was  not 
allowed  to  translate  during  these  tests.  The  specimen  floats 
freely  at  points  "G”  and  "F”  as  point  "F”  deflects.  The 
variable  distance  between  the  load  point  "E"  and  the  actual 
load  point  on  the  specimen,  "  F  " ,  Is  accounted  for  in  the 


results.  The  test  jig  Is  constructed  almost  entirely  of  alu 
m 1 num . 

4.3. 1.3  Test  Procedure 

The  procedure  to  measure  the  strain  on  the  Inside  layer 
Is  as  foil ows : 

First  rings  are  cut  from  the  center  "test  section"  por¬ 
tion  of  the  cylinder.  Half-rings  of  width  3/4"  to  1-3/4" 
were  tested  In  the  Jig  and  could  not  Induce  de 1  am  I nat 1  on . 
Next,  rings  were  cut  from  the  end  "grip  sect 'on"  portions  of 
the  cylinder.  Half-ring  specimens  from  this  thicker  ring 
(see  Table  4.1)  tested  well  In  the  Jig,  delaminating  In  sev¬ 
eral  specimens  in  the  vicinity  of  point  I  in  Figure  4.7  and 
continuing  toward  the  grip  end. 


END 


Figure  4.7  --  Gage  Placement  on  Specimen. 


Four 
from  one 


specimens  were  cut  from  two  3/4  Inch  rings  taken 
end  of  an  outer  tension  specimen.  The  outer 


specimens  were  wound,  from  Inside  to  outside,  as  shown  In 


Table  4.1.  Property  values  for  the  181  glass  clc  n  were 
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taken  from  Military  Handbook  17  .  The  ring  test  to 

determine  residual  stresses  In  the  specimen  were  conducted  as 
described  In  Chapter  3. 

it  was  desired  to  measure  the  strain  in  the  vicinity  of 
the  failure  so  strain  gages  were  placed  as  shown  in  Figure 
4.7.  The  gages  used  were  the  same  variety  used  in  all  previ¬ 
ous  testing  and  the  method  of  placement  Is  detailed  in  Micro¬ 
measurements  Bu'letins  B-129-2^  and  B- 127-9*^. 


Layer 

Th i ckness 

W i nd  Ang 1 e 

Mater i a  1 

1 

0.00746 

90 

S-2/LRF-205 

2 

0.00746 

20 

AS- 4/LRF- 205 

3 

0.00746 

-20 

If 

4 

0.00765 

20 

S-2/LRF-205 

6 

0.00765 

-20 

*1 

6 

0.00746 

90 

•  f 

7 

0.01 

181 

GLASS  CL./LRF-205 

8 

0.0549 

90 

5-2/LRF-205 

TOTAL 

=  0.110 

ile  4.1  -  - 

Lam i nate 

configuration  of 

bu i 1 t-up  ends  on 

tension  specimens. 

Testing  was  done  on  a  60,000  pound  Baldwin  test  machine 
at  a  load  rate  of  approx i mate  1 y  100  pounds  per  minute.  Main¬ 
taining  a  constant  load  rate  was  nearly  impossible  due  to  the 
variability  of  the  specimens'  ability  to  absorb  energy  as 
deflection  progressed. 


4.3. 1.4  Test  Results 

in  the  following  paragraphs,  each  half-ring  delamination 
test  is  discussed  with  relevant  information  and  results.  The 
four  specimens  tested  were  arbitrarily  labeled  A2,  A3,  B2 , 


and  B3  with  the  A  specimens  being  two  halves  of  an  original 


ring  specimen  and  1 ikewise  for  the  B  specimens. 

Specimen  A2  was  loaded  with  as  constant  a  load  with 
respect  to  time  as  possible.  Considering  the  scale  of  the 
load  dial  on  the  60,000  pound  Baldwin  the  load  vs.  time  graph 
of  Figure  4. 8. a  appears  to  be  fairly  linear.  The  rate  of 
strain  with  respect  to  load  did  not  change  appreciably 
throughout  the  test  indicating  a  constant  absorption  of 
energy  by  the  specimen  as  seen  in  Figure  4.8.b.  Specimen  A2 
failed  in  sudden  delamination  at  approximately  179  pounds 
with  extrapolated  strains  of  0.0215  and  0.0170  for  gages  one 
and  two  respectively.  Following  delamination  the  residual 
load  supported  by  the  specimen  was  65  pounds. 

NOTE:  The  extrapolation  was  necessary  due  to  the  Inabil¬ 

ity  to  take  a  reading  at  the  point  of  failure. 

Specimen  A3's  strain  vs.  time  graph  (Figure  4.9.c)  indi¬ 
cates  a  f a i r I y  constant  deflection  as  opposed  to  a  constant 
load  (Figure  4. 9. a).  This  enables  us  to  witness  how  load 
rate  changes  and,  most  importantly,  how  the  load-strain  curve 
changes.  As  the  load-strain  graph  (Figure  4.9.b)  shows 
there  is  a  definite  change  in  the  rate  at  which  strain  occurs 
for  a  given  load.  This  indicates  that  at  loads  above  approx¬ 
imately  80  pounds  this  specimen  requires  less  energy  to 
deflect.  Specimen  A3  failed  in  sudden  delamination  at 
approx i mat e I y  139  pounds  with  extrapolated  strains  of  0.028 
and  0.024  for  gages  one  and  two  respectively.  Following 
del  am  I  nation  the  residual  load  supported  by  the  specimen  was 
55  pounds. 

Specimen  B2  follows  the  same  pattern  as  specimen  A3.  The 


strain-time  graph  (Figure  4.10.c)  shows  a  nearly  linear 
path  and  the  load-time  graph  (Figure  4. 10. a)  shows  an 
obviously  nonlinear  path.  The  apparent  nonlinear  portion  of 
the  time-strain  curve  is  due  to  the  full  scale  setting  of 
25,000  microstrain  on  the  strain  indicator.  Again,  notice 
the  obvious  change  in  slope  of  the  load- strain  curve  in  Fig¬ 
ure  4. 10. b.  Although  the  change  is  less  sudden  there  is  a 
definite  decrease  in  stiffness  of  the  specimen  as  It  begins 
to  take  smaller  increases  in  load  to  achieve  the  same  deflec¬ 
tions.  Specimen  B2  never  failed  in  del  am (nation  as  it 
reached  a  load  of  190  pounds  and  deflected  to  a  point  of 
being  comp  1 e t e 1 y  flat. 

Specimen  B3  was  tested  completely  by  using  the  crank  on 
the  testing  machine  to  load  the  specimen  at  two  levels  of 
constant  deflection  (Figure  4 .  1  1  . c ) .  This  method  allowed  us 
to  witness  pauses  in  the  load  rate  more  accurately.  As  the 
load-time  graph  of  Figure  4. 11. a  shows,  there  were  pauses  at 
70,  110,  130,  and  150  pounds.  As  the  specimen  reached  each 

of  these  four  loads  the  load  actually  dropped  back  10-15 
pounds  before  continuing  under  constant  deflection.  The 
load-strain  curve  of  Figure  4.il.b  better  shows  this  lag  in 
load  required  to  produce  strain.  This  specimen  reached  165 
pounds  before  failing  in  sudden  de I  am i nat i on .  It  reached 
strains  of  approximately  0.030  and  0.026  respectively  for 
gages  one  and  two.  The  specimen  supported  a  residual  load  of 
65  pounds  following  failure. 

Using  the  analysis  above.  Table  4.2  was  obtained  show¬ 
ing  calculated  radial  stresses  for  each  specimen  not  i nc 1 ud- 


1 ng  residual  stresses.  The  calculated  residual  stresses  In 
these  specimens  were  found  to  average  -28.4  +/-  2.0  psi  in 
the  Interlaminar  region  between  the  graphite  and  glass  heli¬ 
cals. 


Spec i men 


Max  Strain 
at  fa i 1 ure 
0.0215 
0.0280 

0.0300 


Final  Inside 
Rad i us  (in) 
2.476 
3.285 

3.652 


Radial  Stress 
at  Failure  (ps 

630.2 

619.2 

597.0 
MEAN=6 15.5 
CoC=0. 027 


Table  4.2  --  Results  of  half-ring  delamination  test. 

4.3. 1.5  Analysis  of  Results 

Table  4.3  gives  the  midplane  stresses  in  the  hoop  (cir¬ 
cumferential  )  direction  corresponding  to  the  measured  strains 
for  specimen  B3  (also  character i st i c  of  specimens  A2  and  A3). 
These  stresses  include  the  calculated  residual  stresses  det¬ 
ermined  earlier  for  the  thickened  specimens. 


Table  4.3  --  Midlayer  stresses  and  loads  through  the 
thickness  of  specimen  B3  at  the  point  of  failure. 

The  stresses  in  Table  4.3  are  very  close  to  failure 


I 


resulting  In  a  de 1  am f nat i on  of  the  specimen  due  to  transverse 
load  and  Interlaminar  shear. 


To  get  rid  of  the  transverse  load  and  thus  the  high 
shear ,  it  is  tempting  to  come  up  with  a  method  to  apply  a 
pure  moment  to  the  half-ring  specimen.  This  would  result  in 
only  norma  1  radial  stresses  in  the  out-of-plane  de 1  ami nat ion- 
causing  stresses.  However,  returning  to  the  theory  used  to 
reduce  the  measured  strain  to  radial  stresses  and  plotting 
calculated  radial  stress  against  measured  strain  for  the 
entire  possible  deflection  (flattening  out)  of  the  half-ring 
specimen,  we  see  that  the  maximum  possible  radial  stress  for 


this  thickness  and  material  stacking  sequence  is  about  630 
ps i  ( f i gure  4.13). 


This  level  of  radial  stress  is  reached  before  delamina¬ 


tion  occurs,  indicating  that  the  radial  stress  allowable  must 


be  higher  than  630  psi.  Another  means  of  determining  this 


allowable  must  be  investigated. 


Using  the  Chang-Spr i nger  failure  criteria  for  fiber 


reinforced  composite  bends,  the  problem  may  become  clear. 


Figure  4.14  indicates  this  failure  envelope  with  the  solid 


line  being  the  envelope  representing  the  commonly  used  val¬ 


ues£  for  interlaminar  shear  and  normal  allowables  in  epoxy 


resin  based  laminates.  The  circles  represent  the  values  of 


shear  and  radial  stresses  calculated  by  the  half-ring  test. 


The  single  square  symbol  is  a  result  obtained  from  the  coupon 


bench  test  discussed  in  the  next  section. 


COMMON  VALUfS 


rCmup*n 


MmtfRini 


INTERLAMINAR  KM  STRESS  R— THETA  (KSI) 


Figure  4.14  --  Chang-Spr  1  nger  Fa i  Ice  Envelope 
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4.3.2  Coupon  Bench  Test 


4.3.2. 1  Introduction 

The  goal  of  this  test  Is  to  apply  a  direct  tens! le  dela- 
rn  i  nat  1  ng  load  to  a  small  section  of  our  specimen.  Lagace  and 
Weems'7  tested  AS4  graphite/epoxy  In  out-of-plane  tension. 

They  used  a  thick  (12-13  mm)  laid  up  laminate  specimen  necked 
down  In  the  test  section.  The  specimen  was  bonded  to  alumi¬ 
num  grips  which  could  provide  stabilization  of  the  specimen 
and  apply  the  out-of-plane  tensile  load.  The  necking  down  of 
the  test  section  a  I  lowed  them  to  produce  a  higher  stress  In 
the  plane  of  desired  f a i lure  than  in  the  plane  of  adhesive 
bonding.  Lagace  and  Weems'7  were  able  to  produce  good  data 
with  an  average  through- the-th i ckness  strength  of  6.24  ksi. 

A  shortcoming  of  this  method  involves  the  fabrication  of 
the  specimen  and  its  thickness.  The  specimen  laminate  may 
not  have  been  laid  up  with  the  same  standards  as  the  materi¬ 
als  for  which  this  type  of  information  is  sought.  Perhaps  test 
specimens  were  laid  up  with  different  thicknesses,  different, 
pressures  app I i ed ,  different  cure  c yc l es .  or  even  different 
f  abr  f  cat  i  on  techniques  such  as  fit  ament,  winding.  These  vari 
ables  could  change  the  through-the  thickness  strength  of 
app I  i c  ab l e  c  ompo  sites. 

The  goa>  of  the  coupon  bench  test  is  to  develop  this  same 
principle  to  work  for  coupons  Cut  from  f *  lament  wound  compos 
i te  cylinders,  namely,  those  be i ng  test ed  here.  Ip**  spe*  < 
mens  wer  e  on  I  y  about  I  .  i  arid  2  .  H  mm  t  h  I  <•  w  t  f  r  t  he  test 
sect  i  on  per  t  i  on  of  t  tie  ♦  es  t  r  ,  <  r  ,Oe  r  and  toe  end  «  e ,  ♦  *  •  '  C  •'  '  .'■] 


tlons  respect f ve 1 y .  This  small  thickness  would  not  a  I  low  a 


necking  down  In  the  "test  section"  of  the  coupon.  Therefore 
failures  near  the  load  discontinuities  under  tensile  loading 
would  be  possible  and  may  cause  premature  failure.  The  value 
of  this  test  Is  that  It  can  provide  a  pure  radial  stress 
allowable  providing  failure  occurs  away  from  harmful  discon¬ 
tinuities  or  If  the  discontinuity  effects  are  small.  If 
failures  occur  In  the  test  section  of  the  coupon,  then  we 
have  a  specific  allowable.  If  failure  occurs  near  the  ends 
where  end  discontinuities  are  In  effect,  the  test  still  pro¬ 
vides  a  useful  1  lower  boundary  to  the  radial  stress  allow¬ 
able. 


Test . 


4. 3.2.2  Testing  Apparatus  and  Procedures 

Figure  4.15  Indicates  the  loading  grips  that  were 
mach  I  ned  spec  I  f  I  ca  I  I  y  for  I  am  1  nates  of  these  rad  I  I  and  thic¬ 
knesses.  The  half-inch  holes  In  each  grip  provide  both  load¬ 
ing  points  and  modi  I  Ity  of  the  specimen  to  reduce  any  pos¬ 
sible  damag i ng  moments  within  the  specimen.  The  grips  were 
machined  From  a'umtnum  with  the  "square"  grip  having  radius 


on  its  concave  surface  equal  to  the  outer  radius  of  the  half- 
r f ng  specimen.  The  round  grip  was  machined  with  a  radius 
equal  to  the  inner  radius  of  the  specimen.  The  thickness  of 
the  grips  was  cut  and  milled  to  0.723  inches.  Specimens  were 
carefully  cut  with  a  bandsaw  to  avoid  fraying  of  the  ends. 
Final  specimens  ranged  from  0.71  to  1.02  inches  wide.  The 
second  planar  dimension  of  the  coupon  was  that  of  the  grips 
since  they  were  sanded  flush  to  the  surface  of  the  Jig  after 
bond i ng . 

Bonding  was  done  in  a  jig  which  supported  both  grips  to  a 
plate  with  a  set  screw  to  tighten  the  grips  together  through 
the  plane  of  the  coupon.  The  coupon  was  centered  in  the  jig 
by  sight  using  hash  marks  scribed  into  the  surface  of  the 
"square"  grip.  3M  Scotch-Weld  DP- 1 1 0  clear  epoxy  adhesive 
was  used  for  bonding.  The  surfaces  were  prepared  using 
methyl  ethyl  ketone  as  a  cleaning  agent  on  the  aluminum  and 
the  normal  surface  preparation  techniques  for  applying  the 
strain  gages  in  earlier  tests.  The  adhesive  is  fast  setting 
and  flexible  with  a  manufacturer’s  claim  of  2500  ps I  shear 
strength.  No  values  of  normal  strength  were  given.  From  48 
to  72  hours  set  time  is  recommended. 

4. 3.2.3  Test  Results 

Five  tests  were  completed  as  time  allowed.  Four  thick¬ 
ened  end  coupons  and  one  thinner  test  section  coupon  were 
tested.  Table  4.4  shows  the  results  of  these  tests  inclu¬ 
ding  failure  descriptions. 


TEST  #  SPEC.  FAILURE 
TYPE  OESCR. 


END 

END 

END 

END 

MIDDLE 


REA 
1  n  ) 


0.737 

0.690 

0.657 

0.713 

0.513 


STRESS 


FAILURE  DESCRIPTIONS: 

1.  Failed  under  impact  due  to  error.  Adhesive 
fa i I u re . 

2.  Clean  failure  in  last  outer  glass  hoop  ply. 

3.  Clean  failure  between  inner  glass  hoop  and 
graphite  helicals. 

4.  Clean  failure  between  graphite  and  glass 
he  1 i ca 1 s . 


Table  4.4  --  Results  of  Coupon  Del  ami  nation  test. 


The  failures  in  the  end  specimens  may  exhibit  some  load 

discontinuities,  perhaps  explaining  their  lower  failure 

stresses.  The  failure  in  the  one  test-section  specimen  was 

exactly  where  actual  del  ami  nation  occurred  in  the  original 

test  cylinders.  This  failure  stress  of  2200  psi  can  then  be 
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plotted  on  the  Chang-Spr i nger  failure  envelope  of  Figure 
4.15.  These  results  indicate  a  smaller  failure  envelope  than 
Lagace  and  Weems17  reported.  Differences  in  this  value  among 
various  laminates  is  certainly  not  inconceivable.  The  coupon 
tested  from  the  test  section  of  the  cylinder  had  gaps  in  the 
wfnding  pattern  of  the  graphite  helicals  of  approximate  1 y  15% 
as  viewed  in  a  Zeiss  stereo  microscope.  These  gaps  certainly 
effect  the  radial  strength  of  the  material.  The  question 
arises  whether  or  not  one  expects  filament  wound  laminates  to 
have  lower  transverse  strengths  than  laid  up  laminates  due  to 


a  greater  percentage  of  voids,  wind  tension,  or  other  factors 
of  the  fabrication  process.  Questions  like  this  must  be  ans¬ 
wered  before  genera)  allowables  can  be  used  with  confidence. 

Certainly,  more  tests  are  necessary  to  establish  a  data 
base  upon  which  to  draw  conclusions.  This  test,  however, 
does  show  some  promise  as  a  means  of  determining  radial 
stress  allowables  in  composite  cylinders. 


CONCLUSIONS 


1.  The  general  purpose  proprietary  finite  element  pro¬ 
gram  ANSYS  is  a  viable  tool  for  the  three-dimensional  study 
of  laminated  composite  cylinders.  It  is,  however,  costly  and 
time  consuming. 

2.  Linear  coefficients  of  thermal  expansion  for  filament 
wound  laminates  can  be  predicted  using  classical  lamination 
theory . 

3.  Residual  stresses  in  filament  wound  laminates  can  be 
predicted  using  classical  lamination  theory. 

4.  Progress  toward  practical  normal  del  ami  nation 
strength  tests  in  cylindrical  laminates  has  been  achieved. 
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APPENDIX  -  STRAIN  GAGE  INFORMATION 

Gage  Type:  EA-06-250BG- 1 20 

EA  -  General  Purpose  constantan  strain  gage 
( EA  Is  the  series) 

(the  M I cro- Measurement s  EA  gage  Is  like 
the  BLH  Electronics  FAE) 

I  mil  (.03  mm)  polylmlde  backing 
temperature  range  -320  F  to  400  F  for  dynamic 
test  I ng . 

STRAIN  RANGE  =  +/-  5?.  FOR  G.L.  >  1/8” 

E=open- faced  general  purpose,  polyimide  backing 
A=constantan  alloy  in  se 1 f-temperature-compensated  form 
06=the  S-T-C  ( Se I f - Temperat ur e-Compensat i on )  number  is  the 
approximate  thermal  expansion  coefficient  in  PPM/  F  of  the 
structural  material  on  which  the  gage  is  to  be  used.  06  can  be 
for  'A'  al I oys  or  'K'  a! loys. 

250=the  active  gage  length  In  mils 
BG=  specifies  Grid  and  Tab  geometry 
1 20= res i stance  In  Ohms 
gage  length=6.3S  mm 
overal I  length=9.S3  mm 
grid  width=3.1R  mm 
overall  width =3. 18  mm 
mat  r i v  s  f  ?e= I  3 . ?mml  x  S . hmrnW 
See  MM  engineer  i  nq  data  street  for  more  Informat  ion. 
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